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A second black hole candidate in a M31 globular cluster is 
identified with XMM-Newton 
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ABSTRACT 

We use arguments developed in previous work to identify a second black hole candidate 
associated with a M31 globular cluster, Bo 144, on the basis of X-ray spectral and 
timing properties. The 2002 XMM-Newton observation of the associated X-ray source 
(hereafter XBo 144) revealed behaviour that is common to all low-mass X-ray binaries 
(LMXBs) in the low-hard state. Studies have shown that neutron star LMXBs exhibit 
this behaviour at 0.01-1000 keV luminosities ^10% of the Eddington limit (^Edd)- 
However, the unabsorbed 0.3-10 keV of XBo 144 luminosity was ~0.30 ^Edd for a 
1.4 M Q neutron star, and the 0.01-1000 keV luminosity is expected to be >3-7 times 
higher. We therefore identify XBo 144 as a black hole candidate. Furthermore, it is 
the second black hole candidate to be consistent with formation via tidal capture of 
a main sequence donor in a GC; such systems were previously thought non-existent, 
because the donor was thought to be disrupted during the capture process. 
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1 INTRODUCTION 



In lBarnard et al] (|2008T ). we identified the first genuine black 
hole candidate using our X-ray classification m e thod, which 
was developed to the fullest in I Barnard et al.l ((2008) . The 
X-ray source associated with the globular cluster Bo 45 
(hereafter known as XBo 45) exhibited behaviour associated 
with all low mass X-ray binaries (LMXBs) in the low state: 
hard power law emission , and high r.m.s variability (see e.g. 
van der Klidll994lll995h . A recent survey of Galactic neu- 
tron star LMXBs has shown that the low state is observed 
at 0.01-1000 keV luminosities <0.1 Z/Edd f° r systems that 
trace a diagonal transition in colour-colour space from low 
state to high state, and < 0.02 Z/Edd for systems that ex- 
hibit a vertical transition (|Gladstone e t al. 2007); Z/Edd is 
the Eddington luminosity. However, XBo 45 exhibited this 
behaviour at a 0.3-10 keV luminosity ~120 % Eddington 
for a 1.4 Mq neutron star, hence we identified it as a black 
hole candidate. 

XBo 45 is particularly interesting because there has 
been a distinct lack of black holes in globular clusters (GCs), 
and most black hole LMXBs are tho ught to be ejected from 
the cluster. The theoretical work of Kalog era, et al.l (|2004l ) 
showed that one possible formation channel for black hole 
LMXBs is the tidal capture of a main sequence star; indeed, 
it is the most likely c hannel for dense cluste rs, such as those 
with collapsed cores. iKalogera et al.l ([2004) predicted such 
systems to be persistently bright, and inferred from the ab- 
sence of such systems that the prospective donor is disrupted 



during capture. XBo 45 has appeared in all X-ray observa- 
tions covering that region of sky, spanning ~30 years, and 
hence is consistent with theoretical predictions for a system 
formed by tidal capture. 

In this paper, we use the same arguments to identify 
XBo 144 (r2-5 in iKong et al.ll2002l . a = 00:42:59.803 S = 
41:16:06.01) as another black hole candidate. Furthermore, 
XBo 144 is located in a confirmed M31 GC (see the revised 
Bologna Catalo gue V.3.5, March 2008; iGalleti et ail |2004| 
I2005ll2006l [20071 ). 

We present detailed analysis of the 2002 June XMM- 
Newton observation of XBo 144, the deepest of four expo- 
sures made between 2000 and 2002. Results from the other 
observations are consistent with this one. In Section[3]we de- 
scribe the observations and data reduction, then present the 
results of our analysis in Section [3] We discuss our findings 
in Section 2] and draw conclusions in Section [S] 



2 OBSERVATIONS AND DATA ANALYSIS 

We rebuilt the data products for the 2002 June 26 XMM- 
Newton observation of the central region of M31 using ver- 
sion 7.1 of the Science Analysis Software suite (SAS). In 
order to screen for background flares, we followed the rec- 
ommended procedure, finding a small flare near the start 
of the observation; this was removed. We the n synchronised 
the pn and MOS lightcurves as described in lBarnard et all 
|2007l h 
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We extracted pn dStriider et all 120011 ), and MOS 
(MOS1 and MOS2 lTurner et al.ll2001T )~data from from a cir- 
cular source region with radius 20", with a 20" background 
region at a similar off-axis angle, with no point sources and 
on the same CCD chip. A larger background is desirable, 
but was prevented by crowding. Source and background 
lightcurves were constructed from the pn and MOS data 
in 0.3-2.5, 2.5-10 and 0.3-10 keV bands. The lightcurves 
were background subtracted. Source and background spec- 
tra were then created, along with corresponding response 
matrices and ancillary response files. The pn image yielded 
5542 net source photons, while the data from the two MOS 
cameras were combined, yielding 4893 net source photons. 
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3 RESULTS 
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Figure 1. Unfolded 0.3-10 keV pn spectrum of XBo 144 with 
best fit PL emission model. 



3.1 Spectral analysis 

The pn and combined MOS spectra were modeled simul- 
taneously, using XSPEC ver 11.3. Each model consisted of 
an emission model, with line-of-sight absorption and a nor- 
malisation constant that accounts for differences in the pn 
and MOS calibrations, after setting the pn normalisation 
to 1. The parameters of the emission models were forced 
to be the same for the pn and MOS spectra, but free to 
vary. Three emission models were initially chosen: black- 
body (BB), bremsstrahlung (BR), and power law (PL). The 
blackbody and bremsstrahlung models were characterised 
by kT, where T is the temperature and k is the Boltzmann 
constant, while power law emission is characterised by pho- 
ton index, T. Absorption is expressed in terms of Nh, the 
equivalent absorption by neutral hydrogen. 

We first modelled the XBo 144 spectra with single com- 
ponent emission models : BB, BR and PL. We found the pn 
and MOS spectra to be best fitted by the PL model, with 
iV H = 8.1±1.2xl0 20 atom cm" 2 and T = 1.48+0.04; all spec- 
tral models are provided in Table [T] The best fit model is as 
expected for a LMXB in its low state; however, we do not ex- 
pect such behaviour at X-r ay luminosities abov e ~0.1 L^dd 
in the 0.01-1000 keV band (|Barnard et alj2008l ). The 0.3-10 
keV luminosity for our best fit PL model is 5.33±0.03x 10 37 
erg s _1 , which is 0.29 Z/Edd for a 1.4 Mq neutron star, and 
0.20 Lsdd fo r the most massi ve observed neutron star (2.1 
M , see e.g. iNice et alj|2005l ). This is sufficently bright for 
XBo 144 to be a plausible black hole candidate. 

We then modelled the XBo 144 spectra with a two com- 
ponent (BB+PL) emission model. The favoured blackbody 
temperature is 0.0082+0.0016 keV; this tells us that there is 
no blackbody component in the 0.3-10 keV band. In Figs.[T] 
and [2] we compare the PL and BB+PL emission models 
with the 0.3-10 keV pn spectrum, unfolded from the instru- 
mental responses; we see that the blackbody component in 
Fig. [2] is set to have minimal contribution to the emission. 

We also modelled the pn and combined MOS spectra 
separately with a PL emission model, to see if the instru- 
ments are in agreement. The best fit models are also shown 
in Table [T] the pn and MOS spectra are in excellent agree- 
ment. We are therefore confident in our interpretation of the 
X-ray emission from XBo 144. 
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Figure 2. Unfolded 0.3-10 keV pn spectrum of XBo 144 with 
best fit BB+PL emission model. The PL component is repre- 
sented by a solid line, while the BB component is represented by 
a dashed line. 

3.2 Variability studies 

The combined pn+MOS, background-subtracted, 0.3-10 
keV lightcurve of XBol44 is shown in Fig. [3] with 400 s 
binning. The intensity appears to vary in the manner of low 
state LMXBs; however, the uncertainties are large, and the 
best fit line of constant intensity yields xVdof = 630/601, 
i.e. the lightcurve is consistent with being constant. The frac- 
tional r.m.s. variability is 6+4% on time-scales longer than 
100 s; this is consistent with a LMXB in the low state, but 
adds no extra strength to our argum ent. Without the defi - 
nite variability exhibited by XBo 45 l|Barnard et al.| [2008l. 
the case for XBo 144 being a black hole X-ray binary is some- 
what weaker, and rests on its association with the globular 
cluster Bo 144. 

4 DISCUSSION 

The X-ray source associated with the M31 GC Bo 144, XBo 
144, exhibits an emission spectrum that is characteristic of a 
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Table 1. Best fits to XBo 144 pn and MOS spectra with various emission models. Blackbody (BB), bremsstrahlung (BR) and power 
law (PL) emission models were used, all suffering linc-of-sight absorption. For each model, we provide the absorption (iVn / 10 20 atom 
cm -2 ), temperature (kT / keV), photon index (r), and MOS normalisation constant (riMOs)- Numbers in parentheses indicate 90% 
confidence limits. We then provide the x 2 /dof. No uncertainties were calculated for the BB model, as x 2 /dof = 8.5. 
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Figure 3. Combined, background-subtracted pn+MOS 
lightcurve for XBo 144 in the 0.3-10 keV band. The lightcurve is 
binned to 400 s. 



LMXB in the low state: the 0.3-10 keV emi ssion is described 
by a p ure power law, with V = 1.48±0.04. iGladstone et al.l 
(2007) showed that neutron star LMXBs only exhibit this 
behaviour at 0.01-1000 keV luminosities <10% of the Ed- 
dington limit; the 0.3-10 keV luminosity of XBo 144 exceeds 
0.20 LEdd for all known neutron stars (2.1 Mq or less). 

Estimating the 0.01-1000 keV luminosity of XBo 144 is 
not as simple as extrapolating the power law to 1000 keV. 
Low state spectra are thought to be caused by unsaturated 
inverse Compton scattering of cool photons on hot electrons 
in a corona, resulting in a power law for energies lower than 
the elect ron temperature, and a Wie n spectrum for higher 
energies l|Sunvaev fc Titarehuklll98Ch . The coronae in black 
hole X-ray binaries tend to have tem peratures of ~ 100-300 
keV dMcClintock fc Remillardl l200fj ) ; hence we cannot de- 
termine the corona temperature with 0.3-10 keV data. We 
obtained lower limit estimates to the 0.01-1000 keV lumi- 
nosity of XBo 144 by measuring the flux of a power law with 
T=1.48 for energy ranges 0.3-100 keV and and 0.3-300 keV. 
We therefore expect the 0.01-1000 keV luminosity of XBo 
144 to be at least ~3-7 times higher than the observed 0.3- 
10 keV luminosity, i.e. >0.6-1.4 I/Edd for a neutron star with 
mass 2.1 Mq. Hence, we propose XBo 144 as a candidate 
GC black hole LMXB. 

Without the strong time-variability observed in XBo 
45, it is possible that XBo 144 is simply an active galactic 



nucleus (AGN), since Y ~1.4 for typical AGN. Hence, asso- 
ciating the X-ray source with a GC is strong evidence that 
it is an X-ray binary. We therefore calculated the probabil- 
ity of an AGN coinciding with a GC within the observed 
field of view. We employed the X-ray luminosit y functions 
(XLFs) for AGN devised bv lMoretti et al.1 l|2003l ). They cre- 
ated AGN XLFs in 1-2 keV and 2-10 keV bands; there- 
fore we calculated the absorbed 1-2 and 2-10 keV fluxes of 
XBo 144 (1.08xl0~ 13 and 5.69xl0" 13 erg cm -2 s~\ respec- 
tively). From the XLFs of Moretti et al. (2003), we expect 
3.8±0.9xl0~ 7 AGN per square arcsec to exhibit 1-2 keV 
fluxes as bright as XBo 144, and only 8±5xl0~ s AGN per 
squar e arcsec to exhib it 2-10 keV fluxes as bright as XBo 
144. iKong et all i|2002l ) found the Chandra position of XBo 
144 to be 2.2" from the GC position. There are 20 GCs 
within the field of view of the XMM-Newton observation we 
are interested in; hence, the probability for a chance coin- 
cidence between a GC and an AGN as bright as Bo 144 in 
the 2-10 keV band within 2.2" is 2.4±1.5xl0" 5 . 

The 0.3-10 keV lightcurve of XBo 144 is consistent with 
a low state LMXB, but also is consistent with being con- 
stant. We note that XBo 144 is ~5 times fainter than XBo 
45, and may be more severely limited by low photon counts. 
There is a remarkable similarity between the two systems: 
they exhibit very similar emission spectra, they are both 
associated with GCs, and both appear to have been persis- 
tently bright over the ~30 year period of observation. We 
infer that XBo 144 and XBo 45 are the same class of ob- 
ject: black hole LMXBs, likely produced by tidal capture of 
a main sequence star. 

We now consider the influence of metallicity on 
the probability of fi nding black hole LMXBs in GCs. 
iBellazzini et all l| 19951 ) discovered a novel mechanism to pro- 
mote tidal capture at higher metallicities, in addition to the 
previously suspected correlation between metallicity and ini- 
tial mass function (IMF) . They assume that for a fixed clus- 
ter density, the rate of tidal capture depends on the mass and 
radius of the capturing star. They find that higher metallic- 
ity stars have larger masses and radii; therefore, they expect 
the tidal capture rate to increase with metallicity. Further- 
more, a metal-rich star will more easily fill its Roche lobe. 
They conclude that this effect alone could explain the ob- 
served ratio between frequencies of X-ray sources in metal- 
rich and metal-poor clusters, although there is no reason to 
exclude the effects of metallicity on the IMF. 

iFan et all (|2008h have produced the most comprehen- 
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sive survey of metallicities of M31 GCs yet made. They com- 
bine spectroscopically-derived metallicities for 295 GCs and 
GC candidates with colour-derive d metallicities for 209 GCs 
and candidates. I Fan et all (120081 ) found a mean [Fe/H] of 
— 1.29±0.03 for the combined sample, which is remarkably 
similar to th e mean [Fe/H] for the Milky Way (-1.30±0.05 
lHarrisI Il996l , updated in 2003). In particular, [Fe/H] = 
— 0.6±0.2 for Bo 144, making it ~5 times more metal rich 
than the mean for M31 or the Milky Way. Bo 45 also has 
a relatively high metallicity: [Fe/H] = — 1.1±0.3, suggesting 
that it is ~60% more metal rich than the mean, with pos- 
sible metallicity range 80-300% of the M31 mean. The high 
metallicity of Bo 144, and the inferred increase in tidal cap- 
ture rate, support the plausibility of XBo 144 being black 
hole LMXB. 

5 CONCLUSIONS 

We have identified a second black hole candidate in a M31 
GC; suc h a system would m ost likely be formed by tidal 
capture. iKalogera et al.l (120041 ) found tidal capture to be the 
dominant mechanism for LMXB formation in collapsed core 
GCs, or eve n denser clus ters such as M15. A survey of Galac- 
tic clusters (|Harrisll 19961 ) showed -20% of clusters to be core 
collapsed; therefore, we expect no more than 20% of GC X- 
ray sources to harbour black holes. Traditional methods of 
identifying a black hole involve obtaining a mass function 
from the radial velocity curves of the donor, and are quite 
unsuitable for identifying black holes in GC LMXBs because 
the donor is near-impossible to identify. Our X-ray identi- 
fication of black holes is not so hindered, but does identify 
only that subset of black hole LMXBs in the low state at an 
observed luminosity >3xl0 37 erg s _1 . Applying our X-ray 
technique to other known, bright X-ray sources in globular 
clusters may reveal further candidate black hole systems. 
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